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We study an approximate renormalization-group trans- 
formation to analyze the breakup of invariant tori for three 
degrees of freedom Hamiltonian systems. The scheme is im- 
plemented for the spiral mean torus. We find numerically that 
the critical surface is the stable manifold of a critical nonpe- 
riodic attractor. We compute scaling exponents associated 
with this fixed set, and find that they can be expected to be 
universal. 

PACS numbers: 05A5.+h, 64.60.Ak 



The breakup of invariant tori is one of the key mecha- 
nism of the transition to chaos in Hamiltonian dynamics. 
For two dimensional systems and for frequencies like the 
golden mean, it has been observed that, at the transition, 
a sequence of periodic orbits approaches geometrically a 
torus of the given frequency, with a nontrivial scaling be- 
havior This self-similarity has been described in 
terms of a nontrivial fixed point of a renormalization- 
group transformation |^-^. The sequence of periodic 
orbits responsible for the breakup is generated by the 
continued fraction expansion of the frequency. For the 
extension to systems with three degrees of freedom (d.f.) 
involving three incommensurate frequencies, we lack a 
theory that generalizes the continued fractions. Numeri- 
cally, three d.f. Hamiltonian systems (or equivalently four 
dimensional volume-preserving maps) have been studied 
with an extension of Greene's criterion |p|-p2|. The con- 
clusion of these analysis was that there is no geometrical 
accumulation of periodic orbits around the critical torus, 
and thus absence of universality. 

The aim of this Letter is to show that one can still ex- 
pect universal behavior in the breakup of invariant tori, 
described by a critical attractor of a renormalization- 
group transformation. The idea is that all Hamiltoni- 
ans attracted by renormalization to this set will display 
sequences of scaling factors that appear in a different 
order but with a universal statistical distribution. The 
dominant unstable Lyapunov exponent characterizes the 
approach to criticality of the universality class. 
Nonperiodic attractors for renormalization maps have 
been conjectured and observed in statistical mechan- 
ics and dynamical systems [|l|,|l6|. In Ref. 0, 
a strange attractor was found for renormalization of cir- 
cle maps, and in Ref. p^ ] similar evidence was found 
for area-preserving maps in two dimensions, from the 
scaling analysis of periodic orbits. The origin of ran- 
domness in these latter studies is due to the randomness 



of the sequence of continued fraction approximants for 
an ensemble of the considered frequencies. In contrast, 
in the present three d.f. case, the rational approximants 
are a regular sequence, obtained by iteration of a single 
unimodular matrix, which allows us to define a renor- 
malization transformation with a fixed frequency vector. 
The idea is to set up a transformation TZ that maps a 
Hamiltonian H into a rescaled Hamiltonian Ti{H), such 
that irrelevant degrees of freedom are eliminated. The 
transformation TZ (which will be defined precisely below) 
should have roughly the following properties: TZ has an 
attractive integrable fixed point Hq (trivial fixed point) 
that has a smooth invariant torus of a given frequency 
ujq. Every Hamiltonian in its domain of attraction V has 
a smooth invariant torus with frequency vector Wq. The 
aim is to show that there is another fixed set A which 
lies on the boundary 91? (the critical surface) and that is 
attractive for every Hamiltonian on dT>. The numerical 
implementation for two d.f. gives support to this 

picture. This is by no means trivial since the construc- 
tion of the renormalization iteration is based on prop- 
erties that are valid close to Hq. The numerical results 
indicate that the domain of convergence of the iteration 
TZ indeed extends up to the critical surface. A mathemat- 
ical justification of this observation, and the formulation 
of conditions for its validity, are completely open prob- 
lems. 

In this paper, we study the extension of these ideas to 
systems with three frequencies, by analyzing an approxi- 
mate renormalization transformation based on the work 
of Escande and Doveil [p^-p2[ . Its properties are found in 
agreement with the general picture described above. We 
find universal exponents associated with the fixed set A. 
They are universal in the sense that they only depend on 
the frequency vector considered, and not on the chosen 
one-parameter family. These results give new insights 
for the set-up of an exact renormalization scheme in the 
spirit of Refs. |-§. 

The transformation we define acts on the following class 
of Hamiltonians with three d.f., quadratic in the actions 
A = {Ai, A2, A3), and described by three even scalar 
functions of the angles y ~ ((pi, ip2, ^pa)'- 

H{A,^)^Ul + m{^))in-Ay 



[u}o + giv)^] ■ A + f{ip), 



(1) 



where m, g, and / are of zero average. The vector 
is the frequency vector of the considered torus and Jl 
is a vector not parallel to Wo, with norm one: ||r2|| = 



1 



+ 1^2^ + |J^3p)^^^ = 1- This model is the sim- 
plest one involving invariant tori with three frequencies. 
It can be thought as an intermediate case between two 
and three d.f. Its behavior is more complex than in two 
d.f.: For instance, the geometric Aubry-Mather theory 
is, to our knowledge, not available p^ . 
The transformation TZ is defined for a fixed frequency 
vector u^o with three incommensurate components. We 
choose u}q = ((j^,cr, 1), where a « 1.3247 is the spiral 
mean: it satisfies = a + 1 [^,^. Since ujq is Dio- 
phantine, the KAM theorem applies to Hamiltonians (|^) 
(although they are isoenergetically degenerate [^) and 
shows the existence of a torus T with frequency vector 
u}a for a sufficiently small and smooth perturbation con- 
sisting of TO, g and /. The domain of existence of T 
corresponds to a neighborhood of the trivial fixed point 



Ho{A) = {n-Af/2 + u:o-A, 



(2) 



N''-^ui,k> 1} 




for which T is located at A = 0. From some of its proper- 
ties, (T plays a similar role as the golden mean in the two 
d.f. case ||2^. The analogy comes from the fact that one 
can generate rational approximants by iterating a single 
unimodular matrix N. In what follows, we denote reso- 
nance an element of the sequence {f ^ — /vrfc-i 
where Vi = (1,0,0) and 



TV = 



The word resonance refers to the fact that the small de- 
nominators ■ ffc that appear in the perturbation se- 
ries or in the KAM iteration, tend to zero geometrically 
as k increases (wq • i^k — cr^^^ — > as fc — > oo). Our 
hypothesis (which is also the starting point of a gener- 
alization of Greene's criterion in Rcf. ||l^) is that this 
sequence plays a leading role in the breakup of the in- 
variant torus with frequency vector ujf). We build an 
approximate scheme by considering the three main reso- 
nances i^i, V2, and 1/3. The renormalization focuses on 
the next smaller scale represented by the resonances 1^2, 
f3, together with 1/4 = Nv^^ = fi — 1/3. It includes a par- 
tial elimination of the perturbation (the part which can 
be considered nonresonant on the smaller scale, namely 
the mode fi), a shift of the resonances, a rescaling of the 
actions and of the energy, and a translation in the action 
variables. It is, in spirit, close to the type of transforma- 
tion considered in Refs. (7|j|]. 

The approximations involved in this scheme are the two 
main ones used by Escande and Doveil: 

a] A quadratic approximation in the actions: the rescaled 
Hamiltonian TZ{H) is in general, higher than quadratic 
in the actions; in order to remain in the same family of 
Hamiltonians (|l|), we neglect these higher order terms. 

b] A three-resonance approximation: we only keep the 



three main resonances at each iteration of the transfor- 
mation. This is the simplest generalization to three d.f. 
of the Escande-Doveil approach. 

Renormalization transformation. — The approximate 
transformation we define acts on a reduced family of 
Hamiltonians (|l]). We consider the three most relevant 
Fourier modes {ui, i — 1,2, 3} at each step of the trans- 
formation. Hamiltonian (|l]) can be written as 

3 

H{A, if) = Ho{A) + Y,h, {A) cos{v, ■ (^), 

i=l 

where ffp is given by Eq. (||), and hi{A) = m,j.{'D, ■ 
A)y2 + g,^n-A + U^. 

Our transformation combines thus five steps: 
1] A canonical transformation that eliminates the first 
main resonance Ui to the order 0{s). This is performed 
by a Lie transformation Us ■ iv^^) (v'j^Oj S'^^i" 
erated by S{A,ip) = Si{A) sin{ui ■ (f). The Hamilto- 
nian expressed in the new coordinates is given by H' — 
e^p{S)H = H + {S,H} + {S,{S,H}}/2\ + •• •, where 
{ , } is the Poisson bracket between two scalar functions 
of the actions and angles. The generating function S is 
determined by the requirement that the order 0(e) of 
the mode i/i vanishes: {S, Hq} + hi{A) cos(i'i ■ (p) — 0. 
This equation has the solution S{A, ip) — —hi{A) sin(z/'i- 
ip)/Lj{A) ■ vi, where uj{A) = ljq + {Q, ■ A)Q,. This step 
generates arbitrary orders in the action variables. In or- 
der to map the family of Hamiltonians (|l|) into itself, 
we expand H' to quadratic order in the actions, and we 
neglect higher orders. The justification for this approx- 
imation is that as the torus is located at A = for Hq, 
one can expect that for small e, it is close to A = 0. We 
notice that h2{A) and /13(A) are not changed to the or- 
der 0{£^). Furthermore, we neglect all the Fourier modes 
except 0, V2, V3, and V4. We expand the Hamiltonian 
to the order 0{e^). This leads to the expression of H': 



H' = H 



h2 €03(1/2 ■ v) + hi cos(i^3 • (f) 

-{{S,h,})/2 + {S, /13}, 



(3) 



where () denotes the mean value defined as {h){A) — 
Jrj^s h{A, ip) dV/(27r)^ The last term of Eq. (|) contains 
the Fourier mode f 4 of amplitude 

/14(A) ^ [Siui ■ dhs + h^u^ ■ dSi)/2, (4) 

where d denotes the derivative with respect to A. We 
expand to quadratic order in the actions. 
2] A shift of the resonances Vk '-^ i^k-i- a linear trans- 
formation {A,(p) 1-^ {N^^ A, Nip), where N is N trans- 
posed. This step changes the frequency into Nljq = 
a^^Uo (since is an eigenvector of N by construction). 
3] We rcscalc the energy (or equivalently the time) by 
a factor cr, in order to keep the frequency fixed at ljq. 
The vector fl is changed into NQ. We define the im- 
age n' of nhyn' = Nn/\\Nn\\, in order to have n' of 



2 



unit norm. We rewrite the mean-value term of Eq. (g) 
as ({5, hi}) = fi{n ■ A)^ + an- A + const. 
4] In order to keep the image of g with zero mean-value, 
we eliminate a by a translation in the action variables 
A 1-^ A + a (of order O(e^) in the ft direction). The con- 
stant part of the quadratic term of the resulting Hamil- 
tonian is cr||^J7||2(l + ^){n' ■ Af jl. 
5] In order to map this Hamiltonian back into the form 
(|l|), we rescale the actions: H{A,ip) = XH {A/X,<p) , 

withx = a\\Nnf{l + ^i). 

In summary, the transformation is equivalent to a map- 
ping acting on a 11-dimensional space 

I I I tl I I Pi I I fl f~,/\ 

[m^^ , 5^, , , , g^^ , /^^ , m^^ ,g^^J^^,n), 

defined by the following relations 











^cj\\Nn\\g,^^,, 


fl. 
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= /(I + 






^a\\Nn\\h[^'> 
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^Xah^^\ 




n' 


= Nn/\\Nn\\, 





where h]l is the coefficient in {fl ■ Ay of h4 given by Eq. 




Critical surface of the transformation. — The numer- 
ical implementation of this scheme shows that there are 
two main domains separated by a critical surface: one 
where the iteration converges to Hq and the other where 
it diverges to infinity. The renormalization-group picture 
for two-dimensional systems with golden mean frequency 
showed that this surface is the stable manifold of a non- 
trivial fixed point (or nontrivial fixed set related to this 
nontrivial fixed point by symmetries Here, we can- 
not expect any relatively stable nontrivial fixed point: 
this has been highlighted in Rcfs. jl^,^ and can be ex- 
plained by considering the map (^) . This map has no sta- 
ble point (there is only one hyperbolic fixed point which 
corresponds to u)q). The eigenvalues of N are and 
-V^e^*" where a w 27r x 0.1120. The map (|) leads 
asymptotically to a rotation of angle a. As a is close 
to 27r/9, we expect the results to oscillate approximately 
with period 9 as it has been observed in Rcf. ||l^ . Figure 
1 shows the scaling factor A^+g as a function of Xk after 
k iterations on the critical surface. The points near the 
diagonal correspond to approximate period 9 behavior. 
There are however strong deviations from this behavior. 
Figure 2 shows the statistical distribution of the scaling 
factors. 

The iterations on the critical surface converge to a non- 
periodic bounded set A. Figure 3 shows the projection of 



A on the plane {g„-^,m„-^). This set has a codimension 1 
stable manifold, i. e., one expansive direction transverse 
to the critical surface. This set plays, for the system we 
consider, the same role as the nontrivial fixed point of 
the renormalization-group transformation for quadratic 
irrational frequencies in two d.f. Hamiltonian systems. 
In particular, its existence implies universality for one- 
parameter families crossing the critical surface. We de- 
fine exponents that characterize the universality class as- 
sociated with the spiral mean. The mean-rescaling is de- 
fined by A = lim„_,oo (0^=1 ^k)^^^ : where Xk is the value 
of the rescaling after k iterations on the critical surface. 
We also calculate the largest Lyapunov exponent k. The 
result we found is that these limits do not depend on the 
point on the critical surface where we start the iteration 
nor on the initial choice of ft. The coefficients k and A 
depend only on u)q. Numerically, we find n « 0.5870 and 
A « 2.6640. 

In conclusion, the numerical results indicate that for 
the spiral mean, the critical surface of the approximate 
renormalization transformation is the stable manifold of 
a critical set instead of a nontrivial fixed point. This fea- 
ture depends strongly on the characteristics of the eigen- 
values of N, and consequently we do not exclude a pri- 
ori that some specific choices of the frequency vector (or 
equivalently of N) could lead to a nontrivial fixed point 
of a renormalization-group transformation. 
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FIG. 1. Afc+g as a function of Afc. 
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FIG. 2. Distribution of the values of the rescalings A of the 
critical attractor A. 
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FIG. 3. Projection on the plane {g^^ , rriui ) of the critical 
attractor A 
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